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Abstract 

Polygonal chains are fundamental objects in many applications like pattern recognition and 
protein structure alignment. A well-known measure to characterize the similarity of two polyg- 
onal chains is the famous Frechet distance. In this paper, for the first time, we consider the 
Voronoi diagram of polygonal chains in d-dimcnsion {d = 2, 3) under the discrete Frechet dis- 
tance. Given n polygonal chains C in d-dimcnsion (d = 2,3), each with at most k vertices, 
we prove fundamental properties of such a Voronoi diagram VDf{C). Our main results are 
summarized as follows. 

• The combinatorial complexity of VDp{C) is at least r2(nL-5-J); jn fact, even a slice of 
VDf{C) could contain a Loo-metric Voronoi diagram in k dimensions. 

• The combinatorial complexity of VDf{C) is at most 0{Tf^^^'^), for d = 2, 3. 

• Even if each polygonal chain degenerates into a point in three dimensions (3D), e.g., when 
we are given a set P of n co-planar points in 3D, the corresponding three-dimensional 
Voronoi diagram under the discrete Frechet distance, VDf{P), has a combinatorial com- 
plexity of il{n^). 

Keywords: Voronoi diagram, Frechet distance, discrete Frechet distance, combinatorial complex- 
ity, protein structure alignment 
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1 Introduction 



Frechet distance was first defined by Maurice Frechet in 1906 [8j. While known as a famous 
distance measure in the field of mathematics (more specificaUy, abstract spaces), it was Alt and 
Godau who first applied it in measuring the similarity of polygonal curves in 1992 [1]. In general, 
Frechet distance between 2D polygonal chains (polylines) can be computed in polynomial time 
[U [2] , even under translation or rotation (though the running time is much higher) [3l |T6] . While 
computing (approximating) Frechet distance for surfaces is NP-hard [9], it is polynomially solvable 
for restricted surfaces [5]- 

In 1994, Eiter and Mannila defined the discrete Frechet distance between two polygonal chains 
A and B and it turns out that this simplified distance is always realized by two vertices in A and 
B [7J. They also showed that with dynamic programming the discrete Frechet distance between 
them can be computed in 0(|yl||i?|) time. (In fact, the dynamic programming algorithm works 
even when the vertices of the chains are in any fixed dimension.) In [lOj, Indyk defined a similar 
discrete Frechet distance in some metric space and showed how to compute approximate nearest 
neighbors using that distance. 

Recently, Jiang, Xu and Zhu applied the discrete Frechet distance in aligning the backbones of 
proteins (which is called the protein structure- structure alignment problem) [11] . In fact, in this 
application the discrete Frechet distance makes more sense as the backbone of a protein is simply 
a polygonal chain in 3D, with each vertex being the alpha-carbon atom of a residue. So if the 
(continuous) Frechet distance is realized by an alpha-carbon atom and some other point which does 
not represent an atom, it is not meaningful biologically. Jiang, et al. showed that given two planar 
polygonal chains the minimum discrete Frechet distance between them, under both translation and 
rotation, can be computed in polynomial time. They also applied some ideas therein to design 
an efficient heuristic for the original protein structure-structure alignment problem in 3D and the 
empirical results showed that their alignment is more accurate compared with previously known 
solutions. 

On the other hand, a lot is still unknown regarding discrete/continuous Frechet distance. For 
instance, each Frechet distance is a true distance measure, yet nobody has ever studied the Voronoi 
diagram under such an important distance measure. This problem is fundamental, it has potential 
applications, e.g., in protein structure alignment, especially with the ever increasing computational 
power. Imagine that we have some polylines Ai, ^3--- in space. If we can construct the Voronoi 
diagram for ^1,^2, A3,... in space, then given a new polyline B we can easily compute all the 
approximate alignment of B with AiS. The movement of B defines a subspace (each point in the 
subspace represents a copy of B) and if we sample this subspace evenly then all we need to do is 
to locate all these sample points in the Voronoi diagram for AiS. 

Unfortunately, nothing is known about Voronoi diagram under the discrete/continuous Frechet 
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distance, even for the simplest case when each chain degenerates into a point in 3D. In this paper, 
we will present the first set of such results by proving some fundamental properties for both the 
general case and some special case. We believe that these results will be essential for us to design 
efficient algorithms for computing/approximating Voronoi diagram under the Frechet distance. In 
this paper, all of our results are under the discrete Frechet distance unless otherwise specified. 

The paper is organized as follows. In Section 2, we introduce some basic definitions regarding 
Frechet distance and review some known results. In Section 3, we show the combinatorial upper 
and lower bounds for the Voronoi diagram of a set of 2D polylines C, VDp{C). In Section 4, we 
sketch how to generalize the results in Section 3 to a set of 3D polylines. In Section 5, we show the 
quadratic lower bound for a special case where the n polylines degenerate into n co-planar points 
in 3D. In Section 6, we conclude the paper with several open problems. 

2 Preliminaries 

Given two polygonal chains A, B with \ A\ = k and |-B| = I vertices respectively, we aim at measuring 
the similarity of A and B (possibly under translation and rotation) such that their distance is 
minimized under certain measure. Among the various distance measures, the Hausdorff distance 
is known to be better suited for matching two point sets than for matching two polygonal chains; 
the (continuous) Frechet distance is a superior measure for matching two polygonal chains, but it 
is not quite easy to compute p]. 

Let X be the Euclidean plane M'^; let d{a, b) denote the Euclidean distance between two points 
a,b G X. The (continuous) Frechet distance between two parametric curves f : [0,1] ^ X and 
g:[0,l]^X is 

^rif.g) = inf max d{f{a{s)),g{(3{s))), 

a,/3se[0,l] 

where a and (3 range over all continuous non-decreasing real functions with a(0) = /3(0) = and 
a(l)= /?(!) = 10. 

Imagine that a person and a dog walk along two different paths while connected by a leash; 
they always move forward, though at different paces. The minimum possible length of the leash 
is the Frechet distance between the two paths. To compute the Frechet distance between two 
polygonal curves A and B (in the Euclidean plane) of |^| and \B\ vertices, respectively, Alt and 
Godau [1] presented an Od^dj log^(|j4| |i?|)) time algorithm. Later this bound was reduced to 
0(|^||S|log(|^||i?|)) time m- 

We now define the discrete Frechet distance following [7]. 

Definition 2.1 Given a polygonal chain (polyline) in the plane P = < pi, . . . ,pk > of k vertices, 
a m-walk along P partitions the path into m disjoint non-empty subchains {Vi\i=i,,m such that 

^This definition holds in any fixed dimensions. In Section 4, we will cover the 3D case. 
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Vi= < Pki.i+1, ■■■ ,Pki > and = ko < ki < ■ ■ ■ < km = k. 

Given two planar polylines A = < ai, . . . ,0^ > and B = < 61, . . . , 6; >, a paired walk along 
A and B is a m-walk {Ai\i=i,,m along A and a m-walk {^i}i=i..m along B for some m, such that, 
for 1 < i < m, \Ai\ = 1 or \Bi\ = 1 (that is, Ai or Bi contains exactly one vertex). The cost of a 
paired walk W = {{Ai, Bi)} along two paths A and B is 

d}p {A, B) = vasiyi max d{a,h). 

i (a,fe)GAxBi 

The discrete Frechet distance between two polylines A and B is 

dF{A,B) = mmd^(A,5). 
w 

The paired walk that achieves the discrete Frechet distance between two paths A and B is also called 
the Frechet alignment of A and B. 

Consider the scenario in which the person walks along A and the dog along B. Intuitively, the 
definition of the paired walk is based on three cases: 

1. > = 1: the person stays and the dog moves forward; 

2. \Ai\ > \Bi\ = 1: the person moves forward and the dog stays; 

3. = \Bi\ = 1: both the person and the dog move forward. 




Fig. 1. The relationship between discrete and continuous Frechet distances. 

Eiter and Mannila presented a simple dynamic programming algorithm to compute dF{A,B) 
in 0(|^||i?|) = 0{kl) time [7]. Recently, Jiang et al. showed that the minimum discrete Frechet 
distance between A and B under translation can be computed in 0{k^l^\og{k + /)) time, and 
under both translation and rotation it can be computed in 0{kH^\og{k + I)) time They 
are significantly faster than the corresponding bounds for the continuous Frechet distance. For 
continuous Frechet distance, under translation, the current fastest algorithm for computing the 
minimum (continuous) Frechet distance between A,B takes 0{{kl)^{k + ^)^log(A; + I)) time [3]; 
under both translation and rotation, the bound is 0{{k + 1^^ \og{k + /)) time (T6j . 
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We comment that while the discrete Frechet distance could be arbitrarily larger than the corre- 
sponding continuous Prechet distance (e.g., in Fig. 1 (I), they are d(a2, ^2) and d{a2, o) respectively), 
by adding sample points on the polylines, one can easily obtain a close approximation of the con- 
tinuous Frechet distance using the discrete Frechet distance (e.g., one can use d{a2, h) in Fig. 1 (II) 
to approximate (i(a2,o)). This fact was also pointed out by Indyk in [lOj. Moreover, the discrete 
Frechet distance is a more natural measure for matching the geometric shapes of biological se- 
quences such as proteins. As we mentioned in the introduction, in such an application, continuous 
Frechet does not make much sense to biologists. 

In the remaining part of this paper, for the first time, we investigate the Voronoi diagram of 
a set of polygonal chains (polylines) in d-dimension, for d = 2,3. While Voronoi diagram is a 
central structure in geometric computing and has been widely studied [13], it still attracts a lot of 
attention recently [H [12] (including a new annual International Symposium on Voronoi Diagrams in 
Science and Engineering). We hope that our work will facilitate the understanding of the continuous 
and discrete Frechet distance and further enable their applications in various areas, like pattern 
recognition and computational biology. 



3 The Combinatorial Complexity of VDf{C) in 2D 

Let = < ai, a2, . . . , Ofc > and i?; = < 61, 62, • • • , > be two polygonal chains in the plane where 
Oj = {x{ai),y{ai)),bj = {x{bj),y{bj)) and k,l >1. We first have the following lemma, which is easy 
to prove. 

Lemma 3.1 Let A2 = < ai,a2 > cL^d B2 = < ^1,^2 > be two line segments in the plane, then 

dF{A2,B2) = max{d{ai,bi),d{a2,b2)). 

For general polylines, we can generalize the above lemma as follows. 

Lemma 3.2 Let Ak = < ai, 02, . . . , > and Bi = < bi,b2, . . ■ ,bi > be two polygonal chains in 
the plane where k,l > 1. The discrete Frechet distance between A^ and Bi can be computed as 



dFiAk,Bi) = < 



max{d{ai,bi),i = 1,2, . . . ,k} if 1 = 1, 

ma.x{d{ai, bj), j = 1,2, ... ,1} ifk = l, (1) 

^max{d{ak,bi),mm{dF{Ak-i, Bi^i),dF{Ak, Bi^i),dFiAk-i, Bi))} if k,l > 1. 



Proof. As can be seen from Section 2, the discrete Frechet distance can be computed using dynamic 
programming. We assume that the man walks along A^ and the dog walks along Bi. If / = 1 then 
the dog stays at bi and dF{Ak, Bi) = max{d(aj, bi),i = 1,2, ... , A;}. Ilk = 1 then the man stays at 
ai and dF{Ak,Bi) = max{d(ai, bj),j = 1, 2, . . . , 
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Suppose that both k and I are greater than 1. If the man and dog both move in the last step 
then 

dp{Ak,Bi) = ma,x{d{ak,bi),dF{Ak-i, Bi^i)). (2) 
If only the man moves in the last step then 

dpiAk, Bi) = max(d(afc, h), dpiAk^i, Bi)). (3) 
If only the dog moves in the last step then 

dpiAk, Bi) = max(d(afc, h), dpiAk, Bi^i)). (4) 

The smallest value of ([2]), ([3]) and (HI) is 

max{d(afc, bi),mm{dF{Ak-i, Bi_i),dF{Ak, Bi_i), dpiAk-i, Bi))}. 

The lemma follows. □ 
Based on the above lemma, we try to investigate the combinatorial complexity of VDf{C), the 
Voronoi diagram of a set C of n planar polylines each with at most k vertices. Following [6l [15], a 
Voronoi diagram is a minimization of distance functions to the sites (in this case the polylines in 
C). We first quickly review a result on the upper bound of lower envelopes in high dimensions by 
Sharir |15j . 

Let E = {(Ti, . . . , an} be a collection of n (d— 1) -dimensional algebraic surface patches in d-space. 
Let A(T,) be the arrangement of S. The result in holds upon the following three conditions. 

(i) Each cTj is monotone in the xiX2 ■ ■ ■ -direction (i.e. any line parallel to Xf^-axis intersects 
CTj in at most one point). Moreover, each is a portion of a {d — l)-dimensional algebraic surface 
of constant maximum degree b. 

(ii) The projection of cjj in x^-direction onto the hyperplane = is a semi-algebraic set 
defined in terms of a constant number of {d — l)-variate polynomials of constant maximum degree. 

(iii) The surface patches cjj are in general position meaning that the coefficients of the polyno- 
mials defining surfaces and their boundaries are algebraically independent over rationals. 

Theorem 3.1 115^ Let be a collection of n {d — \)- dimensional algebraic surface patches in d- 
space, which satisfy the above conditions (i),(ii), and (iii). Then the number of vertices of AiT,) 
that lie at the lower envelope (i.e., level one) is 0{n'^^^^''), for any e > 0. 

We now try to show a general upper bound on the combinatorial complexity of VDf{C). 

Polyline-point correspondence. Let ci, C2, . . . , be the sequence of vertices of a polygonal chain 
C and let (a;(cj), y(cj)) be the coordinates of vertex Cj,i = 1,... ,A;. Using this notation we view 
every polygonal chain C with k vertices in the plane as a point in M?^ and vice versa. 
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Lemma 3.3 Let B G M he a polygonal chain of I vertices bi,...,bi in the plane where hi = 
{x{bi),y{bi)). Let f : M."^^ ^ M 6e the distance function defined as 

fiC) = dF{C,B), 

where C = (ci, . . . , c^) G M?'^ is a polygonal chain with k vertices and Ci = {x{ci),y{ci)),i = 1, . . . ,k. 
The space M."^^ can he partitioned into at most {kl)l semi-algehraic sets Ri,R2, R3, ■ ■ ■ such that the 
function f{C) with domain restricted to any Ri is algehraic. Thus, the function f{C) satisfies the 
above condition (i) and (ii). 

Proof. We consider (2) • (2) manifolds in defined as 

{x[ci) - x{bj)f + {y{ci) - y{hj)f = [x{ci>) - x{bj,)f + (^(qO - y{hf )f 

for every four integer i, i',j,j' such that 1 < i < i' < k and I < j < j' < I- They partition M^'^' into 
at most (kiy. semi-algebraic sets Ri, R2, R3, ■ ■ ■ corresponding to the order of distances between Cj 
and hj for all 1 < i < k and 1 < j < /. 

Equation ([1]) in Lemma 13.21 implies that the function /(C) restricted to a domain Rm satisfies 

fiC) = ^{x{c.) - x{h,)f + {y{c.) - y{b,)f 

for some pair The lemma follows. □ 
We now prove the following theorem regarding the combinatorial upper bound for VDf{C). 

Theorem 3.2 Let C be a collection of n polygonal chains Ci, . . . ,C„ each with at most k vertices 
in the plane. The combinatorial complexity of the Voronoi diagram \T>p{C) is 0{n^^^'^), for any 
e > 0. 

Proof. Let cn, Cj2, . . . , Cjfc be the sequence of vertices of the chain Cj, i = 1, . . . ,n and let {x{cij),y{cij)) 
be the coordinates of the vertex Cij,j = 1, . . . ,k. Using this notation we view every polygonal chain 
C with k vertices in the plane as a point in M."^^. The Voronoi diagram VDp{C) can be viewed as a 
diagram in R^*^. It is well-known that a Voronoi diagram can be interpreted as a minimization of 
distance functions to the sites [T5] . 

Let fi : R^'^ — > R, i = 1, . . . , n, be the distance function defined as 

f,{C) = dF{C,Ci), 

where C = (ci, . . . , Ck) G R^'^ is a polygonal chain with k vertices and q = (x(cj), y(ci)), i = 1, . . . ,k. 
By Lemma 13.31 the function fi{C) satisfies the conditions (i) and (ii). And condition (iii) is easily 
satisfied using an argument similar to [15J. The Voronoi diagram VDf{C) corresponds to the lower 
envelope in the arrangement of the surfaces z = fi{C) in R^*^+^. By Theorem 13. II the combinatorial 
complexity of VDpiC) is 0{n^''+^). □ 
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Combinatorial lower bound for VDp{C). We now present a general lower bound for VDp{C). 
In fact we show a much stronger result that even a slice of VDf{C) could contain a Loo-Voronoi 
diagram in k dimensions, whose combinatorial complexity is J7(nL~5~J). 

Schaudt and Drysdale proved that a Loo-Voronoi diagram in k dimensions has combinatorial 
complexity of Q{n^~^^) [14]. Let S = {pi,p2, ■ ■ ■ ,Pn} be a set of n points in M.'^ such that the Loo 
Voronoi diagram of S has complexity of r2(n'-~^-l). Let M > be a real number such that the 
hypercube [—M, M]^ contains S and all the Voronoi vertices of the Loo-Voronoi diagram of S. We 
consider a A:-dimensional flat F of M^'^ defined as F = {(ai, M, 02, 2M, . . . , a^, kM) \ ai, . . . ,ak G K} 
and a projection vr : F — > M'^ defined as 7r(6) = (61, 63, ... , b2k~i), for b = (61, b2, 63, b2k~i, ^2fc)- 

Let C = {Ci, C2, C„}, each Ci being a planar polygonal chain with k vertices. Let Cj = 
< Cii,Ci2,...,Cik > and Cim = {x{cim),y{cim)), for m = l,2,...,k. We set Qm = {pim,mM), for 
l<i<n, l<m<A;. Clearly, every Cj G F. With Cj we associate a point C- = 7r(Cj) in M.^. We 
show that the intersection of F and a FDir(C) has complexity of Q{n^~^^). 

Consider a point T ^ F such that T' = vr(T) is a Loo-Voronoi vertex of S in M^. Then 
T' G [— M, M]'^. At this point, the question is: what is the discrete Prechet distance between T 
and a chain Q? Note that d^{T, d) < M if and only if W = Wq where Wq = {(tm, Qm) | = 
l,...,k}. Therefore dF{T,Ci) = d^°{T,Ci) = max{|j;(ti) — x{cii)\,\x{t2) — x{ci2)\, ■■■,\x{tj) — 
x{cij)\, ...,\x{tk) — x{cik)\}- This is exactly the Loo-distance between T' and or dpiT^Ci) = 
doo(T' ,Cl). Then the slice of VDf(C) contains the Loo-Voronoi diagram of S" in A; dimensions. We 
thus have the following theorem. 

Theorem 3.3 The combinatorial complexity 0/ VDi?(C) for a set C of n planar polygonal chains 
with k vertices is Q(n^~i~^); in fact even a k-dimensional slice ofYDpiC) can have a combinatorial 
complexity of Q,{n^^^ ). 

We comment that this lower bound is probably not tight and some significantly different method 
is needed to improve either the lower or upper bounds (or both). Nevertheless, the above theorems 
show the first non-trivial bounds for the size of VDf{C). 

4 The Combinatorial Complexity of VDf{C) in 3D 

For protein-related applications, the input are polygonal chains in 3D. It turns out that Lemmas 
13. II and 13.21 hold for polygonal chains in M^. Similar to Lemma 13.31 we can prove 

Lemma 4.1 Let B G R'^' be a polygonal chain of I vertices bi, . . . ,bi in M'^, where 
bi = {x{bi),y{bi), z{bi)) . Let f : M^'^ ^ M 6e the distance function defined as 

fiC) = dF{C,B), 
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where C = (ci, . . . , Ck) G M^*^ is a polygonal chain with k vertices and Ci = {x{ci),y{ci), z{ci)),i = 
l,...,k. The space M'^'^ can be partitioned into at most {kl)\ semi- algebraic sets Ri, R2, R3, ■ ■ ■ 
such that the function f{C) with domain restricted to any Ri is algebraic. Thus, the function f{C) 
satisfies the above condition (i) and (ii). 

The main difference is that the manifolds in M?^ are defined as 

{x{ci)-x{bj)f+{y{ci)-y{bj)f+{z{ci)-z{bj)f = {x{c,,)-x{bj,)f+{y{ci,)-y{bj,)f+{z{ci,)-z{bj,)f. 

Using Lemma l4.ll we can prove similar to Theorem 13.21 the following bound. 

Theorem 4.1 Let C be a collection of n polygonal chains Ci, . . . ,C„ each with at most k vertices 
in . The combinatorial complexity of the Voronoi diagram \Df{C) is 0{n^^^'^), for any e > 0. 

We comment that the lower bound in Theorem 3.3 is in fact a special case for VDpiC) in M^. 

5 Lower Bound for a Simple Degenerate Case 

In this section, we present a lower bound for a simple degenerate case. Under this degenerate 
case, the geometric bisectors are easy to visualize. Nevertheless, we will see that Voronoi diagram 
under the discrete (and continuous) Frechet distance is much more complex compared with the 
corresponding Euclidean case. In fact, for the same construction the former Voronoi diagram has a 
combinatorial complexity of Q.{n^) while the latter Voronoi diagram (under the Euclidean metric) 
is only of size 0(n). 

First, let all the polygonal chains be 2D line segments, i.e., C = {Ci\Ci = < ai,bi >,ai = 
{aii,ai2),bi = {bii,bi2), for 1 < i < n}. Then we set 5j2 = for all i. As a matter of fact, C is 
equivalent to a set P of n points in 3D, i.e., P = {{an, aj2, foji)|l < « < n}. 

We now further construct points in P in two equal parts, i.e., P = Pi U Qi, with Pi = {pi = 
(0,0, i)|l < i < n/2} and Qi = {(j, 0, 0)|1 < j < n/2} (without loss of generality, assume that n is 
even) . 

Let us look at the bisector between pi and pi+i. Be reminded that pi and Pi+i are in fact 
line segments; i.e., pi = < (0,0), (i,0) > and pi+i = < (0,0), (i + 1,0) >. Let s = {x,y,z,0) 
be a point on the bisector of pi and Pi+i in 3D. Then the discrete Frechet distance dpipi^s) is 
dp{pi,s) = maxjlz — i\, yxP^y^}. Similarly, dF{pi+i,s) = max{|z — i — 1|, y^a;^ + y^}. 

So, what does this bisector look like? When 1^ — z| and |z — i — 1| are bigger than \/ x"^ + y^, 
then the bisector is a plane z = when ^Jx^^^i-y^ is bigger the bisector are two paraboloids 

[z — i)^ = + and {z — i — 1)^ = + y^. The intersection of the plane z = ^^^"^ and these two 
paraboloids is, of course, a circle + y^ = (|)^ (on the plane z = ^^). See Fig. 2 (part I). 



8 




X 



X 



lune 



part I 



part II 



Fig. 2. Bisectors for pi,pi+i and for qj,qj+i. 



Next let us look at the bisector between qj and ^j+i- Again, be reminded that qj and qj+i 
are in fact line segments; i.e., qj = < (j, 0),(0,0) > and qj+i = < (j + 1,0), (0,0) >. For a 
point s = {x,y,z) on the bisector of qj and qj+i in 3D, we have the discrete Frechet distance 
dF{qj,s) = max{|2:|, {x - j f + y2} and dpiqj+i^s) = max{|2:|, {x - j - 1)^ + y^}. 

In this case, the bisector is even more interesting (and complex): When \z\ is bigger, the bisector 
contains all the points satisfying > (x — j)^ + and > {x — j — l)"^ + and this is not even 
a manifold. Geometrically, this is the intersection of the two paraboloids z^ > {x — j)^ + and 
> [x — j — 1)^ + and on each vertical slice z = c, for some constant c, this is a lune determined 
by two disks {x — j)'^ + < and {x — j — 1)^ + < c^. Of course, when \z\ is smaller, the 
bisector is exactly the plane 

X = See Fig. 2 (part II). 

So how do we determine the combinatorial complexity of VDp(Pi U Qi)? In this case, the 
Voronoi vertex is not a geometric point. We prove below that we can take the first intersection 
of the bisector oi pi,pi+i and the bisector of qj,qj+i, which is a geometric point, and then show 
that this intersection point is farther to some points in Pi — {pi,Pi-\-i} and Qi — {qj,qj+i} than 
to pi^pij^i^qj^qjj^i and it is never closer to any point in Pi — and Qi — {qj^qjj^i} than 

to pi,pi^i,qj,qj^i. Then this point contributes an 1^(1) cost to the combinatorial complexity of 
VDf{Pi U Qi). Finally we show that there are 0,(n) such pi,pi^i pairs and Q{n) such qj,qj^i 
pairs. And this will conclude the proof that VDp{P) has a combinatorial complexity of r2(n^), 
or, put it in another way, given n co-planar points in 3D the Voronoi diagram under the discrete 
Frechet distance has a combinatorial complexity of Q{'n?). This is significantly different from the 
corresponding Euclidean Voronoi diagram of n co-planar points, has a combinatorial complexity of 



Now let us finish the technical details. The first intersection of the bisector of pj,pj+i and the 
bisector of gj,gj+i is determined by = (j; — j — l)^ + y^ and {z — i)'^ = x^ + y^. Therefore, we have 



e(n) [T3|. 
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the intersection point Uj = (^^, — — j- , 5). Clearly, if two pairs and are different 

then tij 7^ tiij'. We now show that Uj is never closer to points in Pi — {pi,Pi+i} and Qi — {qj,qj+i} 
than to pi,Pi-^.i,qj,qj+i (and Uj is farther to some points in Pi — {pi,pi^i} and Qi — {qj, qj+i} than 
to pi,pi-^.i,qj,qj+i). This can be done by showing dF{tij,Ph) > | for ^ > z and dF{tij,Ph) = | for 
h < i; and dpiUj, qm) > | for m > j + 1 and dpitij^qk) = | for m < j + 1. We have 

dF{tij,Ph) = maxjy (^— )^ + 1 , I2 - 

which is dpitijjPh) = max{|, || — h\} > |, when /i > i and dF{tij,Ph) = max{|, || — h\} = |, when 
h < i. Similarly, 

c^Fl^ij, gm) = max{y(^ m)^ + ,|-|}, 

which is dpiUj, qm) = max{ ' ^ || > |^ when m > j + 1 and 

dpitijiQin) = max{ ^ -4m(j+i)+4m ^ || ^ |^ when m < j ' + 1. Be reminded that if two pairs 

and are different then t^j 7^ tj/j/. Consequently, each bisector of Pi,Pi+i and each bisector of 

qj,qj+i determines a distinct Voronoi point (on a Voronoi vertex) which is never closer to points 

in Pi — {pi,Pi+i} and Qi — {qj,qj+i} than to pi,pi+i,qj,qj+i and is farther to some points in 

Pi — {pi,Pi+i} and Qi — {qj,qj+i} than to pi,pi+i,qj,qj+i. We thus have the following theorem. 

Theorem 5.1 The combinatorial complexity of YDp{P) for a set P of n degenerate planar line 
segments, which is a set of n co-planar points in 3D, is O(ri^). 

Since the continuous and the discrete Prechet distance are identical for line segments, we also 
have the following corollary. 

Corollary 5.1 The combinatorial complexity of the Voronoi diagram, for a. set P of n degenerate 
planar line segments (which is a set of n co-planar points in 3D) under the continuous Frechet 
distance measure is r2(n^). 

6 Concluding Remarks 

In this paper, for the first time, wc study the Voronoi diagram of polylines in 2D and 3D under 
the discrete Frechet distance. Wc show combinatorial upper and lower bounds for such a Voronoi 
diagram. At this point, to make the diagram useful, we need to present efficient (approximation) 
algorithms to construct such a Voronoi diagram for decent k (say A: = 10 ~ 20) so that one can 
first use a {k — l)-link chain to approximate a general input polyline. Although the running time 
might still look too high (due to the 0{k) exponent in the running time), the good news is that 
in many applications like protein structural alignment, n is not very large. Another open problem, 
obviously, is to improve the combinatorial bounds shown in this paper. 
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